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f*^ Abstract 

The present work deals with the problem of a semi-infinite crack steadily propagating in 
an elastic body subject to plane-strain shear loading. It is assumed that the mechanical 
response of the body is governed by the theory of couple-stress elasticity including also 

^ \ micro-rotational inertial effects. This theory introduces characteristic material lengths in 

~Y~, order to describe the pertinent scale effects that emerge from the underlying microstruc- 

ture and has proved to be very effective for modeling complex microstructured materials. 
It is assumed that the crack propagates at a constant sub-Rayleigh speed. An exact full 
field solution is then obtained based on integral transforms and the Wiener-Hopf tech- 
nique. Numerical results are presented illustrating the dependence of the stress intensity 

>~ factor and the energy release rate upon the propagation velocity and the characteristic 

material lengths in couple-stress elasticity. The present analysis confirms and extends pre- 
vious results within the context of couple- stress elasticity concerning stationary cracks 

i— 1 by including inertial and micro-inertial effects. 
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1 Introduction 

It is well known that classical continuum theories possess no intrinsic length scale and thus 
fail to predict the scale effects observed experimentally in problems with geometric lengths 
comparable to the lengths of the material microstructure (e.g. Fleck and Hutchinson, 1997). 
This discrepancy between the classical theoretical predictions and experimental results is found 
more pronounced for materials with a coarse-grain structure. In fact, the macroscopical be- 
havior of most microstructured materials, like ceramics, composites, cellular materials, foams, 
masonry, bone tissues, glassy and semi-crystalline polymers, is strongly influenced by the 
microstructural characteristic lengths, especially in the presence of large stress (or strain) 
gradients (Maranganti and Sharma, 2007). 

On the other hand, generalized continuum theories intend to capture the effects of mi- 
crostructure by enriching the classical continuum with additional material characteristic length 
scales, and, thus, extending the range of applicability of the 'continuum' concept in an effort to 
bridge the gap between classical continuum theories and atomic-lattice theories. A thorough 
review of generalized continuum theories can be found in Maugin (2010). 

One of the most effective generalized continuum theories has proved to be the theory 
of couple-stress elasticity, also known as Cosserat theory with constrained rotations. This 
theory is the simplest gradient theory in which couple-stresses make their appearance. In 
particular, the couple-stress theory assumes an augmented form of the Euler-Cauchy principle 
with a non-vanishing couple traction, and a strain-energy density that depends upon both 
the strain and the gradient of rotation. Such assumptions are appropriate for materials with 
granular structure, where the interaction between adjacent elements may introduce internal 
moments. In this way, characteristic material lengths may appear representing the material 
microstructure. The presence of these material lengths implies that the couple-stress theory 
encompasses the analytical possibility of size effects, which are absent in the classical theory. 
The fundamental concepts of the couple-stress theory were first introduced by Cauchy (1851), 
Voigt (1887) and the Cosserat brothers (1909), but the subject was generalized and reached 
maturity only in the 1960s through the works of Toupin (1962), Mindlin and Tiersten (1962), 
and Koiter (1964). 

Early applications of couple-stress elasticity dealt mainly with stress-concentration prob- 
lems concerning holes and inclusions. In recent years, the couple-stress theory and related 
gradient theories attracted a renewed and growing interest in dealing with problems of com- 
plex microstructured materials. This is due to the inability of the classical theory to predict 
the observed scale effects, and also due to the increasing demands for manufacturing devices 
at very small scales. This approach and related extensions have been recently employed suc- 
cessfully to model size effects in microstructured materials in, among other areas, fracture (see 
e.g. Huang et al. 1997; Chen et al. 1998; Huang et al., 1999; Shi et al, 2000; Georgiadis, 2003; 
Radi and Gei, 2004; Gourgiotis and Georgiadis 2007; Radi, 2008; Aravas and Giannakopoulos, 
2009; Gourgiotis and Georgiadis, 2011; Piccolroaz et al., 2012; Sciarra and Vidoli, 2012b), 
plasticity (see e.g. Fleck et al., 1994; Gao et al., 1999; Hwang et al., 2002; Dal Corso and 
Willis, 2011), and wave propagation (see e.g. Vardoulakis and Georgiadis, 1997; Georgiadis 
and Velgaki, 2003; Georgiadis et al., 2004; Engelbrecht et al., 2005; Polyzos and Fotiadis, 2012; 
Gourgiotis et al., 2013). 

For materials with microstructure, the characteristic material length mentioned before may 
be on the same order as the length of the microstructure. For instance, Chen et al. (1998) 
developed a continuum model for cellular materials and found that the continuum description 
of these materials obey a gradient elasticity theory of the couple-stress type. In the latter 
study, the intrinsic material length was naturally identified with the cell size. Mora and 
Waas (2000) performing experiments in honeycomb materials estimated that the value of the 



characteristic material length £ in couple-stress elasticity is Wd ~ 15d, where d is the diameter 
of the cell of the honeycomb. Also, Chang et al. (2003) associated the microstructural material 
constants of the couple-stress theory with the particle size and the inter-particle stiffness in 
a granular material. In addition, couple-stress theory was successfully utilized in the past to 
model some materials with microstructure like foams (Lakes, 1983) and porous solids (Lakes, 
1993). Generally, the couple-stress theory is intended to model situations where a material 
with microstructure is deformed in very small volumes, such as in the immediate vicinity of 
crack tips, notches, small holes and inclusions, and in micrometer indentations. A recent 
study by Bigoni and Drugan (2007) provides an interesting account of the determination of 
couple-stress moduli via homogenization of heterogeneous materials. 

Regarding now solutions closely related to our problem, we note that due to the com- 
plexity of the equations of couple-stress theories, very few dynamic crack problems have been 
considered in the literature. In particular, Itou (1972) studied the plane-strain time-harmonic 
problem of a mode I finite-length crack in a Cosserat medium with no micro-inertia employing 
the method of dual integral equations. In addition, Itou (1981) evaluated numerically the 
stress intensity factor for the problem of a propagating Yoffe crack in the context of standard 
couple-stress elasticity, analyzing the influence of the crack tip speed on the asymptotics of 
the stress field. Han et al. (1990) investigated the dynamic propagation of a finite-length 
crack under mode-I loading in a micropolar elastic solid by numerically solving a system of 
singular integral equations. They provided solutions for dynamic stress intensity and cou- 
ple stress intensity factors by using the obtained values of the strengths of the square root 
singularities in macro-rotation and the gradient of micro-rotation at the crack tips. Later, 
Georgiadis (2003) using an analytic function method formulated the time-harmonic problem 
of semi-infinite crack under anti-plane strain conditions in a material exhibiting gradient ef- 
fects. Recently, Mishuris et al. (2012) using the Wiener-Hopf technique solved the problem 
of a semi-infinite mode III crack steadily propagating in an elastic solid using the theory of 
couple-stress elasticity including micro- rotational effects. They showed that the process zone 
near the crack tip is strongly influenced by the microstructural parameters, such as the char- 
acteristic material lengths, the micro-inertia, and the maximum crack speed. Finally, Itou 
(2013) solved problem of two collinear finite-length cracks subjected to time-harmonic stress 
waves impinging normal to the crack faces by using the dynamic couple-stress theory with no 
micro-inertial effects. The Schmidt method was utilized to satisfy the boundary conditions 
along the crack faces. He presented asymptotic results for the stress and the couple-stress 
intensity factors and examined the influence on the solution of the ratio of the characteristic 
material length to the pertinent geometrical lengths of the problem. 

In the present work, the dynamic theory of couple-stress elasticity is employed to deal with 
the problem of a mode II semi-infinite crack propagating steadily at a constant sub-Rayleigh 
speed. In this way, the static analysis presented in Gourgiotis et al. (2012) is extended to the 
case of steady-state propagation in order to study the effects of micro-inertia and crack-tip 
speed on the stress and deformation fields, as well as the variation of the fracture toughness 
due to the presence of the microstructure. A micro-inertia term is included in our study 
because previous experience with couple-stress analyses of surface waves and anti-plane crack 
problems showed that this term is important, especially at high frequencies (Georgiadis and 
Velgaki, 2003; Mishuris et al., 2012). In fact, including this term in the present problem gives 
dispersion curves that mostly resemble with the ones obtained by atomic-lattice considerations. 
It is worth noting that this is the first study in the literature concerning a crack propagation 
problem under plane-strain conditions in the context of couple-stress elasticity including also 
micro-inertial effects. 

The paper is organized as follows: In Section 2 the fully dynamical version of couple-stress 
elasticity theory with micro-inertia is presented. The micro-rotational inertia term is included 



in our formulation by considering the appropriate expression for the kinetic-energy density of 
the material particle. In the next Section, the basic equations in plane-strain are provided and 
the influence of couple-stresses on the propagation of Rayleigh surface waves is studied. Unlike 
the conventional elastic theories, this reasonably simple gradient theory can indeed predict 
the dispersive character of Rayleigh waves in a medium with microstructure. The steady-state 
mode II crack propagation problem is formulated in Section 4. The analytical full-field solution 
is then addressed in Section 5 based on the Fourier transform and the Wiener-Hopf technique. 
In Section 6, we present numerical results for the stresses ahead of the crack-tip and the 
crack-face displacements. In addition, a closed-form expression for the stress intensity factor 
is provided. The dependence of these quantities upon the crack speed and the characteristic 
material lengths is also discussed in detail. Finally, we evaluate the dynamic energy release 
rate (J-integral) for the steady-state propagating mode II crack. To this purpose, a new 
extended definition of the dynamic J-integral is given in the context of couple-stress elasticity. 
It is shown that the J-integral can be determined through the use of distribution theory using 
only asymptotic results and is strongly influenced by the microstructure of the material. 

2 Fundamentals of dynamic couple-stress elasticity 

In this Section, we present the basic elastodynamic equations of couple-stress elasticity. Our 
derivation of basic results relies on the momentum balance laws. Interesting presentations of 
the theory can also be found in the works by Mindlin and Tiersten (1962), Koiter (1964), 
and Muki and Sternberg (1965). However, the latter formulations do not include inertia and 
micro-inertia effects since they are of quasi-static character. 

As mentioned before, couple-stress elasticity assumes that: (i) each material particle has 
three degrees of freedom, (ii) an augmented form of the Euler-Cauchy principle with a non- 
vanishing couple traction prevails, and (iii) the strain-energy density depends upon both strain 
and the gradient of rotation. In addition, the kinetic-energy density T, within a geometrically 
linear theory, takes the following form (Nowacki, 1986) 

T = T^Uqilq + -UJqUg, (1) 

where p is the mass density, / is the micro-rotational inertia, u q is the displacement vector, 
LOq = l/2e qpk d p u k is the rotation vector, d p () = d()/dx p , the superposed dot denotes time 
derivative, and the Latin indices span the range (1,2,3) (indicial notation and summation 
convention is used throughout). The second term in the RHS of Eq.ffTl), involving the spin, 
represents the micro-rotational inertia of the continuum. This term, which is not encountered 
within classical continuum mechanics, reflects the more detailed description of motion in the 
present theory. 

For a body with bounding surface S and volume V, the balance laws for the linear and 
angular momentum read 

i T&> dS+ [ F q dV= [ P Uq dV, (2) 

Js Jv Jv 

/ (e q pkX p Tf> + M (n) ) dS + (e qpk XpF k + C q ) dV = {e qpk x p u k + Iu q ) dV, (3) 

Js ' Jv Jv 

where e qpk is the Levi-Civita alternating symbol, T q is the surface force per unit area, M q 
is the surface moment per unit area, F q is the body force per unit volume, C q is the body 
moment per unit volume, and x p designate the components of the position vector. 



Next, pertinent force-stress and couple-stress tensors are introduced by considering the 
equilibrium of the elementary material tetrahedron and enforcing (2) and (3), respectively. 
The force stress tensor a pq and the couple-stress tensor fi pq (both asymmetric) are then defined 

by 

T^=a pq n p , M^=n m n p , (4) 

where n p are the direction cosines of the outward unit vector n, which is normal to the surface. 
In addition, just like the third Newton's law T^ n ' = — T' _n ' is proved to hold by considering 
the equilibrium of a material 'slice', it can also be proved that M> n ) = — M'~ n '. The couple- 
stresses n pq are expressed in dimensions of [forcelllength]' 1 . Further, a pq can be decomposed 
into its symmetric and anti-symmetric components as follows 

&pq — Tpq + O pq , (O) 

with T P q = T qp and a pq = —a qp , whereas it is advantageous to decompose fi pq into its deviatoric 
and spherical parts in the following manner 

fi pq = m pq + -5 pq [J,kk, (6) 

where p pq = m pq , p pq = (l/3)S p qfikk) an< l $pq 1S the Kronecker delta. Now, with the above 
definitions and the help of the Green-Gauss theorem, one may obtain the stress equations of 
motion. In particular, Eq. (|2j leads to the following force equation 

dpcr pq + F q = pu q , (7) 

which, by virtue of (pi), becomes 

d p T pq + d p a pq + F q = pu q . (8) 

Further, Eq. (I3J) in conjunction with (jTl) leads to the following moment equation 

d p p pq + e qp ka p k + C q = ICbq, (9) 

which can also be written as 

^pqkdipik + ot pq + -e pq kCk = -e pq ktik- (10) 

Combining Eqs. (pb, (fsT) , (10) and by taking into account that VxV^' s ' = 0, yields the single 
equation 

dpT pq + -e qpk d p dimik + -e qpk d p Ck + F q = pu q + -e qp kd p LJ k - (11) 

which constitutes the final equation of motion. 

For the kinematical description of the continuum, the following quantities are defined in 
the framework of the geometrically linear theory 

£ pq = 2^p u q + d q u p ), uj pq = -(d p u q - d q u p ), (12) 

OJq = —€ q pk(J p Uk, K pq = O p U) q , (13) 

where e pq is the strain tensor, uj pq is the rotation tensor, and n pq is the curvature tensor (i.e. 
the gradient of rotation or the curl of the strain) expressed in dimensions of [length] -1 . Notice 
also that (13)2 can alternatively be written as 

K pq = ^ e qikd p diUk = e q ikdi£ pk . (14) 



Equation (14) expresses compatibility for curvature and strain fields. The compatibility equa- 
tions for the strain components are the usual Saint Venant's compatibility equations. Further, 
the identity dkK pq = dkd p UJ q = d p Kk q defines the compatibility equations for the curvature 
components. We notice also that K pp = and, therefore, the curvature tensor has only eight 
independent components. 

The traction boundary conditions, at any point on a smooth boundary or section, consist 
of the following three reduced force-tractions and two tangential couple-tractions (Mindlin and 
Tiersten, 1962; Koiter, 1964) 

P q n) = a pq n P ~ 2 e iP kn P dkm (nn)i R q n) = m pq n p - m {nn) n q , (15) 

where mi nn ) = n p n q rn pq is the normal component of the deviatoric couple-stress tensor m pq . 
We remark that in the case in which edges appear along the boundary, an additional boundary 
condition should be imposed. Indeed, as Koiter (1964) pointed out, a force (per unit length) 
tangential to the edge should be specified according to the relation: Q = 2l m (nn)]]> where |[| 
denotes the jump of the enclosed quantity through the edge. This tangential line load along the 
edge is the counterpart of the concentrated normal force which may be specified at the corner 
of the edge of a Kirchhoff plate or shell. Moreover, it is worth noticing that at first sight, it 
might seem plausible that the surface tractions (i.e. the force-traction and the couple-traction) 
can be prescribed arbitrarily on the external surface of the body through relations Q , which 
stem from the equilibrium of the material tetrahedron. However, in this case, the resulting 
number of six traction boundary conditions (three force-tractions and three couple-tractions) 
would be in contrast with the five geometric boundary conditions that can be imposed (Koiter, 
1964). Indeed, since the rotation vector uj q in couple-stress elasticity is not independent of the 



displacement vector, see Eq. (13li; the normal component of the rotation is fully specified 
by the distribution of tangential displacements over the boundary. Therefore, only the three 
displacement and the two tangential rotation components can be prescribed independently. 
As a consequence, only five surface tractions (i.e. the work conjugates of the above five 
independent kinematical quantities) can be specified at a point of the bounding surface of 



the body, i.e. Eqs. (15). On the contrary, in the Cosserat (micropolar) theory, the traction 
boundary conditions are six since the rotation is fully independent of the displacement vector 
(see e.g. Nowacki, 1986). In the latter case, the tractions can directly be derived from the 
equilibrium of the material tetrahedron, so (J4| are the pertinent traction boundary conditions. 
For a linear and isotropic material behavior, the strain-energy density has the following 
form 

W = W \£pq, Kpq) = T-^E-pp^qq T" l-l£pq£pq ~\~ ^^pq^pq + ^f] KpqKqp, \^J 



where (A, //, rj, rf) are material constants. Then, Eq. (16) leads, through the standard varia- 
tional manner, to the following constitutive equations 



dW dW , 

a (pq) = a = ^Pq £ kk + 2//£pg, flpq = -r = 4,7] K pq + &T] K qp . (17) 

OEpq OKpq 



In view of (17), the moduli (A,//) have the same meaning as the Lame constants of classical 
elasticity theory and are expressed in dimensions of [force][length]~ 2 , whereas the moduli 
(77,77') account for couple-stress effects and are expressed in dimensions of [force]. Further, 
following Mindlin and Tiersten (1962), we assume W to be a positive definite function of its 
arguments, so that 

n' 
3A + 2/x>0, <u>0, 77 >0, -K — < 1. (18) 

V 



Incorporating now the constitutive relations (17) into the equation of motion (11) and using 



the geometric relations (12) and (13), one may obtain the equations of motion in terms of the 
displacement vector 



pV 2 u + (A + ju)V(V • u) - p£ 2 V 2 [V 2 u - V(V • it)] = pu + -V x V 



x il, 



(19) 



where £ = (rj/fi) is a characteristic material length, and the absence of body forces and 
couples is assumed. In the limit £ — >• 0, the Navier-Cauchy equations of classical linear isotropic 



elasticity are recovered from ( 19 ). Indeed, the fact that Eqs. ( 19 ) have an increased order w.r.t. 



their limit case (recall that the Navier-Cauchy equations are PDEs of the second order) and 
the coefficient £ multiplies the higher-order term reveals the singular-perturbation character 
of the couple-stress theory and the emergence of associated boundary-layer effects. 



Next, by taking the divergence and the curl of (19) we obtain the equations governing the 



propagation of dilatation and rotation, respectively 

4V 2 (V • u) = V 



u. 



<4(i - rv^)v 2 (v x u) = (1 - h 2 V 2 )V x u, 



(20) 

(21) 

where c p = [(A + 2p)/p] 1 ' 2 and c s = (p/p) 1 ' 2 are the velocities of the pressure (P) and shear 
(S) waves, respectively, in the classical elasticity theory, and h = (I/Ap) 1 ' 2 is a characteristic 
intrinsic material length associated with the micro-inertia of the material (Mishuris et al., 
2012). An interrelation of the two characteristic microstructural lengths £ and h was given by 
Georgiadis and Velgaki (2003) by comparing the forms of dispersion curves of Rayleigh waves 
in couple-stress theory with micro-rotational inertia with the ones obtained by the discrete 



particle theory (atomic-lattice approach). It is interesting to note that Eq. (20) governing the 
propagation of dilatational waves is the same as in the classical theory. On the other hand, 
unlike the corresponding case of classical elastodynamics, Eq. (21) is of the fourth order. 



This implies that for shear waves, wave signals emitted from a disturbance point propagate 
at different velocities. The last statement can be supported by considering a time-harmonic 
plane wave solution and determining dispersion relations. Specifically, we consider a plane 



wave solution of Eq. (21 ) in the following form 



u = Ad exp i(£ (n ■ x 



wt) 



(22) 



where A denotes the amplitude, (d, n) are unit vectors defining the directions of motion and 
propagation, respectively, x is the position vector, £ is the wavenumber, w is the circular 
frequency of the plane wave, and i 2 = —1. Then, on substituting (22 ) into Eq. (21 ), we obtain 
the following dispersion relation for the shear waves 



w 



2a2\ 



^ 2 (l + ff)(l + ^ 



,2a2\-1 



(23) 



Accordingly, the phase velocity V s of the shear waves in couple stress elasticity takes the 
following form 

Vs = % = p,(l + ^2 )V 2 (1 + ^ 2) -l/2 (24) 



Equation ( 24 ) shows that the propagation velocity of these waves depends on the respective 



wavenumber. Hence, shear waves are dispersive in couple-stress elasticity, while the longitu- 
dinal waves remain non-dispersive as in the classical theory (Toupin, 1962; Graff and Pao, 
1967). 



To investigate further upon the nature of the dispersion relation in couple-stress elasticity, 
we consider the group velocity v = dw/dE, at which the energy propagates in a dispersive 



medium (Achenbach, 1973). In particular, according to (23) and (24), we obtain 

e - h 2 ) Cs e(i + fer 1/2 (i + h 2 er 3/2 . 



v s + 



(25) 



The following three cases are then distinguished: (i) For £ < h, Eq. (25) implies that v s < V s , 



thus the dispersion for shear waves is normal, (ii) For £ > h, we have v s > V s indicating that 
the dispersion is anomalous, (ra) For £ = h or (£, h) — > (i.e. no material microstructure) , the 
shear wave velocity degenerates into the non-dispersive velocity of classical elastodynamics. 



3 Basic equations in plane-strain 

For a body that occupies a domain in the (x, y)-plane under conditions of plane strain, the 
displacement field takes the general form 

u x = u x (x, y) / 0, u y = u y (x, y) / 0, u z = 0. (26) 

Further, Muki and Sternberg (1965) showed that, if the general 3D equations of Section 2 are 



combined with the restrictions in (26) and the normalization ^k = of the couple-stress field 



is adopted, all field quantities are independent of the coordinate z. Accordingly, except for 
lo z = uj and (k xz , k vz ), all other components of the rotation vector and the curvature tensor 
vanish identically in the particular case of plane-strain considered here. In view of the above, 
the following kinematic relations are obtained 



£.T.T 



"l^H 



- uy 



OyUy, 



1 



£xy 



LO 



yO x lly OyV>x)i 



'UX 



<9 X W, 



(d X Uy +8yU X ), 



x yz 



dyUJ. 



(27) 



(28) 



Accordingly, the constitutive equations in (17) furnish the non- vanishing components of the 



symmetric part of stress and the couple-stress, respectively. Vanishing body forces and body 
couples are assumed in what follows. Then, the antisymmetric part of the stresses is found 



tj.r 



from (10). In view of the above, the following expressions are written 

T~xx — \A ~r ^fJ'jE-xx ~r ^yyi Tyy — V^ ' ^l^l^-yy ' ^xxi 
Tfl xz — Q/J,l K xz , fTlyz — 4/it Ky Z , 

a xx = a yy = 0, a yx = -a xy = 2/j,£ 2 X7 2 uj 
and consequently the (asymmetric) stresses are given as 

Tyx t (Xyxi 



>xy 



2fie 



■>-u 



I 



to-. 



(T.i 



'yy 



'yy 



<J 



!)■>' 



a 



.ri) 



r xy + a 



•»'.'/• 



(29) 
(30) 
(31) 

(32) 



The equation of motion in ( 19 ) leads to the following system of coupled PDEs of the fourth 
order for the displacement components (u x ,u y ) 



fiV 2 u x + (A + fi)d x e - [i£ 2 V 2 [V 2 u x - d x e] = pu x + ^d y u, 
pSJ 2 u y + (A + [i)d y e — fi£ 2 V 2 [V 2 u y — d y e] = pu y 



i 



d x ij, 



(33) 

(34) 



where V 2 = d 2 () + d 2 () and e = d x u x + d y u y is the dilatation. Although the above system is 
much more complicated than that in the respective case of classical elastodynamics, uncoupling 
by the use of Lame-type potentials still proves to be successful. The potentials (j)(x, y, t) and 
ip(x,y,t) are defined in terms of the displacement components as 



d(p dtp dcf) dip 

dx dy' y dy dx 

Consequently, the dilation and the rotation vector become 



(35) 



V^ 



UJ 



-vV 



(36) 



Incorporating the above relations into (33) and (34), we derive, after some tedious but straight- 
forward algebra, the following uncoupled PDEs, which constitute our field equations. The 
potential <p is governed by a second order PDE, whereas the potential ip satisfies a PDE of the 
fourth order, i.e. 

0, (37) 



c 2 V 2 < 



(1 



f 2 V 2 ) 



vV 



;i - /i 2 v 2 ) $. 



(38) 



The fact that ( 35 )- ( 38 ) comprise the complete solution of the equations of motion ( 33 ) and (134 



may be proved by following exactly the steps in Sternberg's (1960) form of the proof for the case 
£ = and h = 0. In fact, it is only necessary to replace in Sternberg's proof the standard wave 
operator governing the shear wave propagation: D 2 , = V 2 — c~ 2 d 2 by the modified operator 



in couple-stress elasticity with micro-inertia: 



01 



[1 - e 2 v 2 ) v 2 - c- 2 (i - h 2 v 2 ) d\ 



Accordingly, the stresses and couple-stresses take the following form in terms of the Lame 
potentials 

72. 



O", 



AV Z 



2// 



a 



y->- 



(7 



■r,V 



/i 



/i 



<Jyy = XV (p + 2fl 

d 2 4> d 2 ip 



+ 



d 2 <p 
dx 2 

dy 2 



+ 



&■ 



dxdy dy 2 dx 2 

d 2 4> d 2 ip d 2 ip 
dxdy dy 2 dx 2 

m xz = -2ii£ 2 d x V 2 iP, 



dxdy J 

d 2 ip \ 
dxdy J 



/urv 4 ^ + 7 v 2 ^, 



+ pLtV A Tp--V 2 lp, 



m 



i-r- 



-2ii£ 2 d y V 2 ^. 



(39) 

(40) 
(41) 

(42) 

(43) 
(44) 



Regarding the traction boundary conditions in the plane strain case, we note that these are 



defined through Eqs. ( 15 ) by taking also into account that the normal component of the couple- 
stress mr nn -j is zero. Indeed, since the components (m xx ,myy,m yx ,m X y) of the couple-stress 
tensor vanish identically in the plane-strain case (recall that in this case, all field quantities are 
independent upon the z coordinate), we conclude that m< nn -> = n p n q mpq = 0. This, in turn, 
implies that no edge forces Q = (l/2)[[m( nn J appear at the corners of a boundary or section 
in plane strain. However, we remark that these edge forces should be considered in antiplane 
strain problems, where, in general, mi nn \ ^ 0. In addition, pertinent edge conditions should 
be taken into account in both plane and antiplane strain problems, in the more general theory 
of dipolar gradient elasticity (see e.g. Gourgiotis et al., 2010; Sciarra and Vidoli, 2012a). 
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3.1 Propagation of Rayleigh waves in couple-stress elasticity 



We close this Section with a brief study of the propagation of Rayleigh waves in couple-stress 
elasticity with micro-rotational inertia. A more detailed study can be found in Georgiadis and 
Velgaki (2003). To this purpose, we consider the following two-dimensional time-harmonic 
response in the half-space y > 



4> (x, y, t) = $(y) • exp [i(£x - wt)] , 
■0 (» 5 V, t) = *(y) • exp [i(£x - wt)] . 



Then, substitution of the above equations into (J3TJ) and (38) yields 

d 2 <S> 



dy 2 



m = 0, 



d 4 ^ 



{Pl + llf* ' ° 2 ' 2 



dy 4 v ^ s ls ' dy 2 
where the following definitions have been employed 



+ #<tf* = 0, 



pi = e 



pi = e 



% 



C 2 _l ^ 2 



with 



with 



w 



1 



V2£ 

1 
V 7 ^ . 



A, 



w 



¥\1 1/2 



A s + 1 



w 2 h? 



1/2 



and 



A, 



1 



u; 2 ^ 2 



■2,„2 



+ 



4£ 2 u; 



1/2 



(45) 
(46) 

(47) 
(48) 

(49) 

(50) 
(51) 

(52) 



The criterion for surface waves is that the displacement potentials decay exponentially with 
the distance y from the free surface. It is noted that in the context of classical elasticity these 



waves are non-dispersive at any frequency. From the solution of the ODEs (47) and (48), we 
obtain the displacement potentials ((f), ijj) as 



i> 



4> = A exp [— P p y + i(£,x — wt)] , 
B exp [—/3 s y + i(£x - wt)] + C exp [-j s y + i(£x - wt)] . 



(53) 
(54) 



Since A s > always, the quantities (a s ,T s ) defined in (50) and (51) are real and greater than 



zero. Consequently, it can be readily shown that 7 2 is always greater than zero, whereas /3 2 
and /3 2 have real values greater than zero if only the phase velocity of the Rayleigh waves 
is c p h = w/£ < min{ V s , c p }, where V^ is the shear wave velocity for a couple-stress material 
defined in Eq. (J24]). 

Now, the pertinent boundary conditions for a traction-free boundary defined by the plane 
(x, y = 0) with n = (0, 1) take the following form 



<Jyy(x, 0) = 0, (Ty X (x, 0) = 0, 



in 



IJZ 



(x,0) = 0. 



(55) 



Substitution of expressions (53) and (54) into the boundary conditions (55) results in an eigen- 



value problem for the wave amplitudes (A, B, C). Accordingly, the vanishing of the pertinent 
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determinant provides the dispersion equation for the propagation of Rayleigh waves in a ma- 
terial characterized by couple-stress elasticity. The analytical expression of the determinant 
is very lengthy and is omitted here. A numerical solution was derived using the symbolic 
computer program MATHEMATICA™. 

Figure Ilk shows the variation of the normalized phase velocity c p h/c s in couple-stress 
elasticity versus the normalized wavenumber £^ = £^ °f the Rayleigh wave for a material with 
Poisson's ratio v = 0.3. It is observed that for small wavenumbers (low frequencies) the phase 
velocity of the Rayleigh wave tends to the respective classical value: cr = 0.927c s . This is 
to be expected intuitively since for relatively long wavelengths the wave should not 'see' the 
material microstructure. On the other hand, for large wave numbers (high frequencies) the 
dispersive Rayleigh wave speed depends strongly upon the ratio of the characteristic material 
lengths ho = h/t. In particular, it is observed that as £d increases, the phase velocity reaches 
a plateau attaining a constant value which depends upon ho and the Poisson's ratio v. Indeed, 
as it is illustrated in Figure [lp, the following three cases are distinguished as ^ — > oo: (i) 
ho < 0.5; in this case the limit value of the normalized phase velocity depends only upon the 
Poisson ratio, (ii) 0.5 < /iq < l/v2; in this case c p h/c s depends upon both ho and v. In fact, 



for v > v*(ho) the ratio becomes constant and tends to the value l//io- (hi) ho 

this case we obtain that c p h/c s = 1/ho always, independently of the Poisson's ratio. Moreover, 
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(b) 



Figure 1: (a) Dispersive character of Rayleigh waves in couple-stress elasticity, (b) Limit of the Rayleigh 
velocity in couple-stress elasticity with respect to the Poisson's ratio as £d — >• oo. 

it is worth noting that the curve /iq = c s /cr in Figure [la, differs slightly from the classical 
elasticity solution, which, in turn, implies that in this case the Rayleigh waves are almost non- 
dispersive. The constant /iq depends, according to its definition, only upon the Poisson's ratio; 
in fact for a material with < v < 0.5, we have 0.874 < cr/c s < 0.955, and, consequently, h^ 
ranges from 1.046 < hg < 1.144 (see e.g. Achenbach, 1973). Finally, it should be remarked 
that a propagation velocity V is characterized as sub-Rayleigh, which is our case of interest 
here, when V < c p h- Consequently, when < ho < /iq it suffices that the propagation velocity 
is less than the classical Rayleigh velocity i.e. V < cr, whereas for ho > h^ the propagation is 
sub-Rayleigh provided that V < c s /ho- 
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4 Formulation of the crack problem 

Consider now a semi-infinite crack in a body of infinite extent under plane strain conditions. 
The body is governed by the equations of couple-stress elasticity. The crack propagates with 
constant sub-Rayleigh velocity straight along the x-axis and is subjected to a distribution of 
shear stresses along the crack faces moving with the same velocity. We now introduce the 
standard steady-state assumption for moving sources according to which a steady stress and 
displacement field is created in the medium w.r.t. an observer situated in a frame of reference 
attached to the tip of the crack, if the crack has been moving steadily (with a velocity V, say) 
for a sufficiently long time. In this way, any transients can reasonably be avoided (therefore 
gaining considerable simplification in the analysis). Specifically, upon introducing the Galilean 
transformation 

X = x - Vt, Y = y, (56) 

the time derivative in the fixed Cartesian system becomes: dt = —Vdx, which, in turn, 
implies that the field and the boundary conditions become independent of the time t, and the 
variables (x, t) enter the problem only in the combination x — Vt. Furthermore, in the new 
moving Cartesian coordinate system, partial derivatives w.r.t. t are neglected and the field 
equations (37) and (38) for the Lame potentials <ft(X, Y) and ip (X,Y) can now be written as 



(1 - 77lV) 



8 2 c 



+ 



o 2 



dX 2 dY 2 



0, 



{ } dx 2 + dY 2 l 



l — m h. 



2l - 2 o) 



d 2 (v 2 vO a 2 (v 2 vO 



ox 2 



+ 



dY 2 



0, 



(57) 



(58) 



1 /9 

where c = c s /c p = [(1 — 1v) jl (1 — z/)] ' < 1. Also, m = V/c s and mc = V/c p are the two 
Mach numbers. It is worth noting that m < tjir with ttir = min{l//\o, l//ig}, always, in order 
for the crack to propagate with a sub-Rayleigh speed (see also Section 2). Figure R^l depicts 
the sub-Rayleigh and the super- Ray leigh regimes in couple-stress elasticity. 



Super-Rayleigh regime 



v = 0.5 




Figure 2: Sub-Rayleigh and Super-Rayleigh regimes in the m — ho plane. 
Due to the anti-symmetry with respect to the Y = plane, the problem can be viewed as a 
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half-plane problem in the region Y > under the following boundary conditions 

a yx (X,Y = 0) = -T(X) 

a yy (X,Y = 0) = 
m yz (X,Y = 0) = 
u x (X,Y = 0) = 

where T(AT) is the distribution of shear tractions along the crack faces. 



for 


- oo <X < 0, 


for 


— oo <X < oo, 


for 


— oo <X < oo, 


for 


<X < oo, 



(59) 
(60) 
(61) 
(62) 



5 Full field solution 

An exact solution of the boundary value problem described above will be obtained here through 
the Fourier transform and the Wiener-Hopf technique (Noble, 1958). The direct and inverse 
Fourier transforms are defined as 



f(s,Y) 



+oo 



f(X,Y)e ls * dX, f(X,Y) 



1 

2^ 



f(s,Y)e-^ds, 



(63) 



where L denotes the inversion path within the region of analyticity of the function in the 



complex s-plane. Transforming the field equations (57) and (58) with (63 )i, we obtain the 
following ODEs 

.72.-, 

(1 - m 2 c 2 )s 2 (j)* = 0, (64) 



~dY* 



£ 



2 dV 

dY* 



1 + (2 - m 2 hl) ts 



d 2 ip* 
~dY^ 



+ 



(l - m 2 ) + (1 - m 2 hl) £ 



2f,2\ n2 2 



2 /* 

S ip 



0. 



(65) 



The above equations have the following general solutions that will be required to be bounded 
as Y — > +oo 

4>{s,Y) = A(s)e- aY , (66) 

iP(s,Y) = B(s)e-P ¥ + C(s)e-~ fY , (67) 



where 



a = a(z) 



[(1 



m 2 c 2^ .211/2 



V]- 



1/2 



P = P(z) 



7 = l(z) 



[l+{2-m 2 h 2 )z 2 + X (z)} 
V2£ 

[l+(2-m 2 h 2 )z 2 - X (z)] l/2 



V2£ 



X = X (z) = [1 + 2 (2 - h 2 ) m 2 z 2 + m% 



4, 4 41 V 2 



with z = si being a dimensionless complex variable. 

The complex function x( z ) nas f° ur branch points at ±ib± and ±162, where 



(68) 
(69) 

(70) 
(71) 



h: 



2-h 2 T 2(l-h$) 1/ ^ i 2 



rnhg 



(72) 



In particular, when ho < 1, 6^2 are always real and the branch points are located along the 
imaginary axis. In this case, the branch cuts are chosen to run from ±z6i to ±262 (see Fig A. la 
in the Appendix A). On the other hand, when ho > 1, 61,2 are complex and the branch points 
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are symmetrically located with respect to the real axis at the four quadrants of the complex 
z-plane (Fig. A. lb). Accordingly, z = and z = ±ibo with 



(1 



m 



2\V2 



'1 



m^ 



,l/2 : 



(73) 



are additional branch points of the functions fi{z) or 7(2) depending on the values of the 
parameters m and h$ (the branch cuts of these functions are given in Appendix A). Note, 
that 6o is always real when m < rriR. In any case, the specific introduction of the branch cuts 
secures that the functions (Xj/?>7)> are single- valued and positive along the real axis. 

Moreover, the transformed expressions for the stresses and displacements that enter the 
boundary conditions take the following form 



< x (s,Y) 



-l> 



dY 4 



l + (2 



m 



- 2 hl)fs 2 ) 



d 2 x/j* 

'dY 2 ~ 



1 - (1 - m z h^)ts z )s^* + 2is 



dY 



v*(a,Y) = (\ + 2fi) 



d 2 <f>* 
dY* 



AS <p + 2lflS 



dY' 



m*(s,Y) 



**(& 



dY 



*( V\ d ^ ■ A* 

Ux(s,Y) = — -iscj) 



(74) 

(75) 
(76) 
(77) 



Next, in preparation for formulating a Wiener-Hopf equation, the unilateral Fourier transforms 
of the unknown stress o~ yx (X > 0, Y = 0) ahead of the crack tip, and the unknown crack-face 
displacement u x (X <0, Y = 0) are defined as follows 



E+(s) = I a yx (X,Y = 0)e lsX dX, a yx (X, Y = 0) = -^ f E' 



(s)e~ lsJi ds, (78) 



and 



17" (s) = / u x (X,Y = 0)e isX dX, u x (X,Y = 0) = ±- f U~ (s)e~ 



isX 



ds, (79) 



where the inversion path is considered to lie inside the region of analyticity of each transformed 
function. In particular, we assume the following finiteness conditions to hold: \a yx (X, Y = 
0)| < M • exp(-p s X) for X -> +oo and \u x {X,Y = 0)| < A • exp( Pu X) for X -> -oo, 
where (M,N,p-£,pu) are positive constants. Consequently, the transformed function E + (s) is 
analytic and defined in the upper half-plane Im(s) > — ps, while U~(s) is analytic and defined 
in the lower half-plane Im(s) < p\j. 



Enforcing the boundary conditions (60)- (62), taking into account also ( 75 )-(|77b and (79)i, 



results in the following equations for the unknown functions A(s), B(s) and C(s) 

21 



A{s) 



m 2 s 



IJ- 



B(s) 



C{s) 



\r 



_(2-m 2 )( 7 2 -s 2 ) 
m 2 /3(/3 2 - 7 2 ) 

(2-m 2 )(p 2 -s 2 ) LT _ 
m 2 7(/3 2 — 7 2 ) 



(80) 
(81) 
(82) 
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Consequently, Eqs. (80)-(82), together with (59) and (74) provide the final Wiener- Hopf 
equation of the problem, connecting the two unknown functions E + (s) and U~(s) 



2\2 „2 



E+(.)-r-(*) 



fj,(2 — m ) s 
m? [{stj 



211/2 



■K(s£)U'(s), 



(83) 



where T (s) is the Fourier transform of the loading along the crack-faces. The kernel function 
K(s£) is given as 



K(z) 



where 



02 + [ z 2]l/2 + m 2 

= 6{z) = £ 



id, 



2 P(zMz) 



(84) 



(85) 



r 2i 1/2 ' 

and d = 4 (l — m 2 c 2 ) / (2 — m 2 ) . The problem has now been reduced to the determination 
of the unknown functions S + (s) and U~(s) from the single functional equation (83). 

5.1 Wiener-Hopf factorization 

To proceed further a product-factorization of the kernel K{z) is required. First it is checked 
that K{z) has no zeros and no poles in the complex plane in the sub-Rayleigh regime (m < 
run). This was verified by using the symbolic computer program MATHEMATICA™. 
Moreover, the kernel function exhibits the following asymptotic behavior 



lim K(z) =£(l-d) + 0(z- 2 ), 

\z\— >oo 



lim K(z) 
|*|-M) 



(1 - m 2 ) 



-1/2 



d 



+ 0(z 2 ). 



(86) 
(87) 



This leads us to introduce a modified kernel given as N(z) = K(z)/ {£{1 — d)], which possesses 
the desired asymptotic property hmM-^oo N(z) = 1. Indeed, this new form of the kernel 
facilitates its product splitting by the use of Cauchy's integral theorem (Noble, 1958; Roos, 
1969). The functional equation (83) takes now the form 

Hl(l-d)(2-m 2 ) 2 s 2 



where the function 
upper and lower ha 



E + (s)-T-(s 

1/2 



m? [si] 1 / 2 ■ [st] 1 ! 2 



N(s£)U~(s), 



(88) 



-[ 



21 V2 



is written as a product of two analytic functions in the 



f-plane, respectively (see e.g. Mishuris et al., 2012) 



[ 



.21 1/2 



,1/2 r ,1/2 



In addition, the modified kernel splits up as 

N(z) = N + (z)-N-(z), 
where 



A^ (z) = exp 



'c d C-z 

I f ln[A(C)] 



2vri 



(- u 



c 



d(l 



(89) 

(90) 

(91) 
(92) 
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The use of the Cauchy integral theorem is depicted in Figure pi The functions N + (z) and 
N~(z) are analytic and nonzero in the half-planes Im(z) > —e and Im(z) < e, respectively, 
with e being a real number such that e — > 0. In fact, introducing e facilitates the introduction 
of the branch cuts for N(z) (see e.g. Georgiadis, 2003). The original integration paths Cj 
and C u extend parallel to the real axis in the complex z-plane. According to the Cauchy 
integral theorem, Jordan's lemma, and by taking into account that N(z) has no poles or 
zeros in the finite complex plane, we are allowed to use equivalent integration paths, i.e. the 
{C' d ,C' u ) contours around the branch cuts of N(z) extending along: ±£ < Im(z) < ±60, and 
±60 < Im(z) < ±61 (Fig. p3j). Note that the second branch of N(z) exists only in certain cases, 
depending on the values of (m, ho) (more details are given in the Appendix B). 



ilm(z) 




Figure 3: Contour integration for the factorization of the kernel function N(z). 



This eventually leads to the following forms of the sectionally analytic functions N (z) 
(the analytic derivation is provided in Appendix B). 



AP^z) =exp 



1 



7T 



ib (] 



tan 



ImiV(C) 
ReiV(C) 



d( f ih _ x 

— \- a tan 

C ± z Jib 



'lmJV(C)" 
_ReA^(C)_ 


dC 

(±z_ 



(93) 



with the properties N + (— z) = N~(z) and N^(z) = 1 for \z\ — > 00. The constant a takes the 
values or 1 depending on the branches of the kernel function N(z), and is defined in Appendix 
B. Moreover, it is noted that when z approaches the branch cuts in the z-plane from above or 
below, the integrals defining the functions N + (z) or N~(z) become singular. However, these 
integrals are not singular simultaneously. Thus, we may calculate these exceptional cases 
without resorting to principal values by employing N (z) = N(z)/N^(z), where we choose in 
the denominator the N (z) function that is not defined by a singular integral. 



In view of the above, Eq. (88) can now be rewritten as 



E+(a) [si] 1 / 2 _ id(l - d) (2 - m 2 ) 2 s 2 



N+(s£) 



m 2 [s£]_ 



N-(s£)U-(s) + 



N+(sl) ' 



(94) 



5.2 Solution of the Wiener-Hopf equation 

We now assume the following form of the loading applied on the crack faces 



T(X) = -± e x/L with - 00 < X < 0, 

1j 



(95) 
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where To and L are positive constants having pertinent dimensions. Transforming (95) with 

To 



(63)i we obtain 



T~ s) 



1 + iaL' 



(96) 



so that (|94|) becomes 

f (*) M 

N+(s£) 



S+(s) [si] 1 / 2 _ fd(l ~d){2- m 2 f s 2 



m 2 \s£] 



N + {s£)(l + isL) 



The sum-splitting of the second term in the RHS of ( 97 ) is required to complete the decoupling 
process. This is be obtained by inspection as 



where 



M(z) 



M + (z) 



To [z 



1/2 



N+(z){l + izt- l L) 



M + (z) + M-{z), 



(l + izl- 1 !;) 



,1/2 



[i£/L] 



1/2 



M~{z) 



N+(z) N+(i£/L) 
[it/L] 1 ' 2 



(98) 



(99) 



(100) 



(l + izl- x L)N+(il/Ly 

with M + {z) being an analytic function in the same upper half-plane where N + (z) is defined, 



while M (z) is an analytic function in the half-plane where Im(z) < i£/L. Equations (99) 



and (100), when combined, allow the final re-arrangement of the Wiener-Hopf equation 

[s£tl 2 



N+(s£) 



M+(s£) 



[i£{\ -d)(2- m 2 ) 



2\ 2 „2 



m 2 [s£] 



1/2 



-N-(s£)U-(s) + M-(s£) = E(s£). (101) 



The above functional equation defines the function E{z) only on the real line. In order to 
evaluate this function, it is first necessary to examine the asymptotic behavior of the functions 
E + (s) and U~(s). In particular, guided by the results concerning the modification of stress 
singularities in the presence of couple stresses in plane-strain crack problems subjected to 
a quasi-static shear loading (Huang et al., 1997; Gourgiotis et al., 2011), we assume that 
the stress and crack-face displacement along the crack faces exhibit the following asymptotic 
behavior 

a yx (X,Y = 0) = O(X-^ 2 ) 

u x (X,Y = 0) = O((-X) 1 / 2 ) 



as 



as 



X - 
X 



+0, 
► -0. 



(102) 
(103) 



FT 



•K-l 

FT 



Further, we consider at this point the transformation formula X K -R- i k+1 T(K + 1) • s 

(with k 7^ 0, +1, +2, ...), where T() is the Gamma function with k > —1. The symbol X4- 
means that the quantities on either side of the arrow are connected through the unilateral 
Fourier-transform. Then, employing theorems of the Abel-Tauber type (see e.g. Roos, 1969), 
we obtain the following asymptotic behavior in the transform domain 



£+(s) = 0( S -V2) 

U-(s) = 0(s- 3 / 2 ) 



as 



as 



+oo with Im(s) > 0, 
+oo with Im(s) < 0. 



(104) 
(105) 



In light of the above, and bearing in mind that: N^(s£) — > 1 and M^(s£) — > as \s\ — > +oo, 



we conclude that the first member of (101) is a bounded nonzero analytic function at infinity 



for Im(s) > 0, whereas the second member is a bounded nonzero analytic function at infinity 
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for Im(s) < 0. Then, based on the theorems of analytic continuation, the two members define 
one and the same analytic function E{z) over the entire complex z-plane. Moreover, Liouville's 
theorem leads to the conclusion that E(z) = Eq, where Eq is a constant. 
The transformed shear stress is now given by ( |101 ) as 

[E + M + (s£)} N + (s£) [s£}~ 1/2 . 



£+(*) 



(106) 



The constant Eq can be determined from simple equilibrium considerations. In particular, 
from the requirement that the upper part of the body Y > is in equilibrium in the moving 
framework, we have: f_ a yx (X, Y = 0) dX = 0, which, by taking into account Eqs. (78 )i 
and (96), can be written as — Tq + S + (0) = 0. Further, in view of (99) and (106), we obtain 



the limit 



E+(0)= lim S + (s) = r 1/2 



M-M) 



m 



-1/2 



\-d 



1/2 



1/2 



En 



Tp [i£/L] 
N+(i£/L) 



-1/2 



+ T + O(s 



1/21 



(107) 
Thus, the only possibility left from ( 107) is the vanishing of the term inside the second bracket, 

T [il/ltf 2 
Eo - N+(il/L) ■ (1 ° 8) 

Alternatively, the constant Eq can be evaluated from the conditions at remote regions in the 



physical plane. First, we observe that the exponentially decaying tractions in (95) can be 
replaced by an equivalent concentrated load of intensity: J_ a yx (X,Y = 0) dX = —Tq at 
distance L from the crack-tip. Thus, according to Saint- Venant's principle and in view of the 
classical solution of the steady-state problem of concentrated shear tractions along the crack 
faces (see e.g. Freund, 1990), we anticipate that for X — > +oo the shear stress ahead of the 
crack-tip will behave as in the classical theory i.e. as ~ X~ 3 ' 2 . This, in turn, implies the 
following asymptotic behavior in the transformed domain: S + (s) = 0(s 1 ' 2 ) as |s| — > 0. Using 



this result in conjuction with Eq. (107) provides the same value for the constant Eq. 



The final transformed expressions (valid for all s in the pertinent half-plane of convergence) 
for the stress ahead of the crack-tip and the crack-face displacement then become 



IT (a) 



£+(*) = 

111 


T [i£/L}\ 


li 


N+{i£/L 

2 


) 


n£(l ~ d) 


[2-m 2 ) 2 





+ M + (sl 



1/2 



Tp [i£/L] 
N+(i£/L) 



N + {s£) [s£} + 1/2 , Im(s) > 
[s£] 1 ! 2 



M~(s£) 



s 2 N-{s£)' 



Im(s) < 0. 



The limits of the latter expressions for \s\ — > oo are found to be 



To [i# 2 
'L-N+iit/L) 

T m 2 



.-1/2 



+ 0(s- 



3/2n 



,1/2 



U (S) fi^L(l-d)(2-m 2 ) 2 N+(i£/LY 
which, accordingly, provide the following near-tip field as 
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6 Results 

6.1 Analytical representation of displacements, stresses and couple-stresses 

The shear stress ahead of the crack tip and the crack opening displacement can be obtained 



from (109) and (110), respectively, by employing the inverse Fourier transform according to 



(63)2- For a crack propagating with a sub-Rayleigh speed, the functions in (109) and (110) 



do not have branch points or poles along the real line, consequently, the path of integration L 
coincides with the real line s. In this case we obtain 



a 



!!■'■ 



(X) 



u x {X) 



1 

1 

2^ 



Z + (s)e- lXs ds, X > 0. 



U-(s)e- tAs ds, X<0, 



(115) 



(116) 



Alternatively, taking into account that S + (s) and U~(s) are analytic in the upper and lower 
half- planes respectively, we can deform the original integration paths by making use of Jordan's 
lemma, around the pertinent branch cuts of these functions. In particular, the branch cuts for 
the functions U~(s) and X + (s) can be chosen to be along the negative (Im(s) < 0) and the 
positive (Im(s) > 0) imaginary axis, respectively. 
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Figure 4: Variation of the crack sliding displacement along the upper crack face for three different values: 
ho — {0, ho, 1.1/iq} and various crack speeds. 

The variations of the normalized tangential displacement on the upper crack face Y = +0, 
X < 0, versus normalized distance X/t to the crack tip are shown in Fig. Hi for a material 
with L/i = 10 and a selected range of values of m, ho and v. It is noted that the specific value 
Ljt = 10 was chosen for most of the graphs in this paper, since for small values of this ratio the 
couple-stress effects are more pronounced. In fact, the effect of the ratio L/l on the dynamic 
solution is the same as in the stationary mode II crack studied previously by Gourgiotis et al. 
(2012). Now, as is in classical elasticity, the crack tip profile is blunted in the couple-stress 
solution. Moreover, it observed, that the magnitude of the sliding displacement between the 
crack-faces increases as the propagation speed approaches the pertinent Rayleigh velocity, for 
every set of material parameters. However, the displacement in couple-stress elasticity remains 
always bounded in the sub-Rayleigh regime. This result is in marked contrast with the classical 
elastodynamic solution, where the crack face displacements become infinite in magnitude as 
we approach the Rayleigh speed (Freund, 1990). This is illustrated more clearly in Fig. p] 
where the variation of the crack face displacement at the point X = —5£ is shown with respect 
to the normalized propagation speed for ho = {0, 0.8/ig}. Since, ho < h^, in both cases, the 
Rayleigh velocity takes the value mji = 1//iq = 0.927 for v = 0.3. As noted earlier, when 
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Figure 5: Variation of the crack sliding displacement at the point X = — 5£ in couple-stress elasticity and 
classical elasticity with respect to the normalized crack speed. 

m —> ttir, the displacement in classical elasticity becomes unbounded, whereas the respective 
displacement in the couple-stress theory increases significantly but remains finite. In addition, 
we note that for every propagating velocity m, the crack-face displacements are smaller than 
the respective ones in the classical theory. This finding shows that the crack becomes stiffer 
due to the effects of the microstructure, here brought into play by rotational gradients. This 
trend was also observed for the stationary mode II crack in couple-stress elasticity (Huang et 
al., 1999; Gourgiotis and Georgiadis, 2007), and proves that the microstructure may shield 
the crack tip from fracture. Finally, we remark that the influence of the Poisson's ratio v on 
the solution becomes significant as the speed of the crack increases, and is more pronounced 
as v — >• 0. 
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Figure 6: Variation of the shear stress ahead of the crack-tip in couple-stress elasticity for three different 
values: ho = {0,/io, l.l/io} and various crack speeds. 

Next, the variation of the normalized shear stress ahead of the crack-tip with respect 
to the normalized distance X/l is depicted in Fig. pi for a material with L/l = 10, ho = 
{0, /iq, 1.1/iq} and v = {0,0.5}. It is observed that the effect of the Poisson's ratio becomes 
more significant when h$ — > (i.e. in the case of no micro-inertia). Further, it is remarked 
that as m — > the steady state solution approaches the static solution, irrespectively of the 
value of the ratio ho, confirming thus previous results obtained by Gourgiotis et al. (2012) for 
the stationary crack. On the other hand, as m — > tur, the zone within which the couple-stress 
effects become important increases. Interestingly, when the propagation speed approaches the 
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Rayleigh velocity, the shear stress takes on negative values very close to the crack-tip within 
the range 2£ < X < 6£, as is shown in Figure ^k for a material with ho = 0.8/ig (thr = 0.927) 
and v = 0.3. However, these negative values are very small in magnitude and appear only in 
the case m — > wir and /iq < h^. In Figure [7k the near-tip couple-stress asymptotic solution 



in (113) and the classical elasticity solution are also depicted. It is noticed that the near-tip 



asymptotic field dominates within a zone of 0.5£ to the crack tip, but the couple-stress effects 
are significant within a zone of 2£. Moreover, as we move away from the crack-tip {X > 8£) 
the full field solution approaches gradually the classical elastodynamic solution for every crack 
speed. 
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Figure 7: (a) Variation of the shear stress ahead of the crack-tip in couple-stress elasticity as the crack speed 
approaches the Rayleigh velocity. The near-tip asymptotic field and the classical elasticity solution are also 
shown, (b) Variation of the couple-stress m xz ahead of the crack-tip for a material with v = 0.3 and ho — 0.8/t ( * 
and various crack speeds. 



Finally, Figure [7b shows distribution the couple-stress m xz normalized by To, versus the 
distance X/£ to the crack tip for a material with v = 0.3 and ho = 0.8/ig- It is observed 
that the couple stress decays rapidly as the distance from the crack-tip increases. Also, it 
is interesting to note that m xz is bounded at the crack-tip, as in the stationary crack case 
(Huang et al., 1997; Gourgiotis and Georgiadis, 2007). In particular, the couple-stress m xz 
takes a finite negative value at the crack tip then increases to a bounded positive maximum, 
and diminishes monotonically to zero for X >> £, thus recovering the classical solution of 
linear elasticity. The couple-stress effects are more significant as the crack propagation speed 
approaches the pertinent Rayleigh velocity ran. 



6.2 The dynamic stress intensity factor 



According to Eq.( 113) and the definition Ku = \im.x^+o(2irX) ' a yx (X,Y = 0), the dynamic 



(d) stress intensity factor (SIF) in couple-stress elasticity becomes 



Kfj 



V2T 



L ■ N+(i£/L) 



;ii7) 



In the limit case m 



0, the dynamic SIF (117) degenerates to its static counterpart K s n 

(2012). The details of the 



in couple-stress elasticity derived previously by Gourgiotis et al. 

pertinent limit procedure and an explicit expression for Kfj are given in Appendix C. 
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Furthermore, the ratio of the dynamic SIF in couple-stress elasticity to the respective one 
in the classical theory (cZ.) is 



where KJf = V2T /VZ. 



K d n _ 1 
K dci. N+(t£/L)' 



(118) 
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Figure 8: Variation of the ratio of the dynamic SIF to the static SIF in couple-stress elasticity with respect 
to normalized crack-speed. 



It is noted that in the classical elastodynamic theory the distribution of the shear stress 
ahead of the crack-tip does not depend upon the crack speed V. Consequently, the dynamic 
stress intensity factor is identical to the corresponding static result for any V: Kjf /Kjj L = 1 
(Freund, 1990). This result is in marked contrast with the one obtained here in the context 



couple-stress elasticity, where the dynamic SIF in (117) depends upon the speed of the propa- 



gation through the function N + (i£/L). An analogous result was found also by Itou (1981) for 
a propagating Yoffe crack in couple-stress elasticity with no micro-inertia. This dependence is 
clearly depicted in Fig. p^for a material with L/£ = 10 and v = 0.3 (/iq = 1.078). Indeed, for 
curves 1-3, we have ho < /iq, which, in turn, implies that the normalized crack velocity should 
be m < 1/Hq in order for the crack to propagate with a sub-Rayleigh speed (see also Fig. Y3. 
On the other hand, when ho > h^ (curves 4 and 5), the requirement for the crack-speed to be 
sub-Rayleigh ism< l//io- It is noteworthy that as the crack speed approaches the pertinent 
Rayleigh limit: m — > txir, a significant but finite increase of the ratio Kfj/ Kfj is observed. 
This increase is more pronounced as ho — > /ij), where Kfj = 9A7Kfj (curve 3). 

Figure [9] illustrates the effect of Poisson's ratio on the variation of the the ratio of the SIFs 
/Kfj, for two values of ho. In particular, in Fig. Km we have ho = 0.8 < /iq for all v 
(recall that 1.046 < /iq < 1.144) and, consequently, m < 1//iq for the crack to propagate with 
sub-Rayleigh speed. This requirement is always satisfied for curves 1-3, however, for m = 0.9 
(curve 4), this inequality is satisfied only in the range 0.137 < v < 0.5. In fact, as v — > 0.137 
the propagation speed approaches the Rayleigh speed, resulting to a steep increase of the 
dynamic SIF. On the other hand, for propagation velocities m < 0.5 the ratio varies slowly. 
Further, in Fig. [9J3 we have ho = 1.2 > /iq for all u, and thus, m < l//io- As m — > 1/ho an 
increase in the ratio Kfj /Kfj is noted again, however, the increase is more moderate compared 
to the case ho < h%. 
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Figure 9: Variation of the ratio of the dynamic SIF to the static SIF in couple-stress elasticity with respect 
to the Poisson's ratio v for: (a) ho — 0.8 and (b) ho — 1.2. 



6.3 The dynamic energy release rate 

In this Section, we evaluate the dynamic J-integral (energy release rate) in the context of 
couple-stress elasticity for a steady-state propagating mode II crack. The J-integral in the 
quasistatic case was first established by Atkinson and Leppington (1974) (see also Atkinson 
and Leppington, 1977; Lubarda and Markenscoff, 2000). By following relative concepts from 
Rice (1968a, b) and Freund (1990), the extended definition of the dynamic J-integral in couple- 
stress elasticity takes the following form 

where T is a piecewise smooth simple 2D contour surrounding the crack-tip, V is the crack- 



tip velocity, W is the strain-energy density (16), T is the kinetic energy density including 
micro-rotational inertia terms defined in M, and (Pq ,Rq ) are the reduced force-traction 



and tangential couple-traction defined in (15). It is noted that in the general transient case 



the J-integral is not path-independent (Freund, 1990). However, for the case of steady crack 
growth considered here, it can be readily proved following an analogous procedure as the one 
proposed by Grentzelou and Georgiadis (2008) for the stationary crack problem, that the J- 
integral remains path-independent. Moreover, we note that in the expression for the dynamic 



J-integral in (119) the edge forces are not included since the latter are zero under plane-strain 



conditions (see also Section 3). Utilizing the Galilean tranformation in (56), the J-integral in 
the steady-state case becomes 
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(120) 



dl 



For the evaluation of the J-integral, we consider the rectangular-shaped contour Y (surrounding 



the crack-tip) with vanishing 'height' along the y-direction and with e — > (see Fig. 10). Such 
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a contour was first introduced by Freund (1972) in examining the energy flux into the tip of a 
rapidly extending crack and it is particularly convenient in computing energy quantities in the 
vicinity of crack tips. In fact, this type of contour permits using solely the asymptotic near- 
tip stress and displacement fields. It is noted that upon this choice of contour, the integral 



J r [VF + T]dY in ( 120 ) becomes zero if we allow the 'height' of the rectangle to vanish. In this 
way, the expression for the J-integral becomes 
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-21im 
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dUq 
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(121) 
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Figure 10: Rectangular-shaped contour surrounding the crack-tip. 

Further, we note that due to the anti-symmetry conditions that prevail in the mode II case, 
the normal stress a yy (X, Y = 0) and the couple-stress m yz (X, Y = 0) are zero along the whole 
crack line. In view of the above, the dynamic J-integral for the mode II plane-strain case 
assumes the following form 

ld or [ e \ <v ,^ 3%(^+°) , (v ^ du x (X,+0) 

J = - 2 }™o] | ayy( - ' +0) ' ~dx — + ayx ( ' + 0) dx — 



+ rn yz (X, + 0). d -^^\dX, 
-2 1imJ^X,+0). du *f^ +0) dX. 



(122) 



Now, by using the asymptotic solution (jl 13|) and (|114|), we finally obtain 

-l/2 (X _)-l/2 dX; 
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2T 2 m 2 



lim 



fmL(d - 1) (2 - m 2 Y [N+(i£/L)] z ^+o J_ e 



(X, 



(123) 



where for any real A with the exception of A = —1,-2,-3,..., the following definitions of 
singular distributions (of the bisection type) have been employed 



Xi 



X x , for X > 
0, for X < 



and 



X ; 



0, for X > 

|X|\ for X < 



(124) 



Moreover, the product of distributions inside the integral in (|123[) is obtained through the use 

(125) 



of Fisher's theorem (Fisher, 1971), i.e. the operational relation: 

(A:_) A (A: + )- 1 - A = -7T(5(X)[2sin(7rA)] _1 with A/ -1,-2, -3, 
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with 5(X) being the Dirac delta distribution. Then, in taking into account the properties of 
the Dirac delta distribution, we obtain the final expression for the J-integral 



J" 



T$m 2 



\xL[d - 1) (2 - m 2 ) 2 [N+(i£/L)] d 



(126) 



Based on the above analysis, we were able now to evaluate the dynamic J-integral. It can 
be readily shown that as I — > 0, the energy release rate (ERR) in couple-stress elasticity 



degenerates into its classical dynamic counterpart. Indeed, by taking into account (87) and 



(90) we obtain 



J. 



cl. 



T^m? 



(1-m 2 ) 



(127) 



1/2 



4(l-m 2 Y /z (l 



m c 



fiL R(m) 
'" 2 ) — (2 — m 2 ) is the Rayleigh function (see e.g. 



where R(m) 
Freund, 1990). 

It is noted that as the speed of the crack reaches the classical Rayleigh wave velocity (i.e. 
rriR = l//ig = c r/ c s) the Rayleigh function becomes zero and, consequently, the ERR in (127) 
becomes infinite. This is a common feature of the steady-state problems in classical elasticity 
where the SIF does not depend upon the crack-tip velocity. On the other hand, for £ ^ (i.e. 



when couple-stress effects are taken into account) the denominator in ( 126 ) is never zero in the 
range m £ [0, iur], which, in turn, implies that the ERR in couple-stress elasticity, remains 
always finite even when the crack propagates with the Rayleigh speed. 
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Figure 11: Distribution of the dynamic ERRs in couple-stress elasticity and classical elasticity with respect 
to the normalized crack-tip velocity for three different values: ho — {0, ho, 1.05/t ( *}. 



Figure [TT] depicts the variation of the normalized ERR in couple-stress elasticity and in 
classical elasticity with respect to the normalized crack propagation velocity m for a material 
with L/£ = 10 and v = 0.3. For curves 1 and 2 the Rayleigh velocity limit is mn = 1//iq = 
0.927 with Hq = 1.078, whereas for curve 3 (ho = 1-1/iq) we have: m^ = 1/ho. As the 
propagation speed approaches the pertinent Rayleigh limit, an increase is observed in the 
dynamic ERR in couple-stress elasticity, however, contrary to the classical elasticity case 
(dashed curve), the ERR remains always bounded. The increase becomes more significant 
when ho = hg (curve 2). This is attributed to the fact that in this case the Rayleigh waves 
travel almost non-dispersively in an infinite medium governed by couple-stress elasticity, thus, 
resembling the classical elasticity situation (see Fig. [T] and the relevant discussion in Section 
3). Finally, applying the classical Griffith criterion for the stability of the crack propagation: 
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J d = J d ri where J d r is the critical value depending on the material propreties, we conclude 
that, for J d r constant, the steady state motion studied here is likely to be unstable and that the 
crack will accelerate rapidly up to the pertinent limiting Rayleigh velocity as in the classical 
theory of elasticity (see e.g. Willis, 1971). 
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Figure 12: Variation of the ratio of the dynamic ERR in couple-stress elasticity to the ERR in classical 
elasticity with respect to the normalized crack-tip velocity: (a) for ho < /lo an d (b) for ho > hg. 

The variation of the ratio of the ERR in couple-stress elasticity to the respective one in 
the classical theory versus the normalized velocity is shown in Figures 12 1 and [j~2p for the 
cases: ho < h$ and ho > h\ 



, respectively. In particular, in the case ho < /ig, the ratio of the 

In fact, 
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as the 



ERRs is a monotonically decreasing function of the crack speed (Fig. 
crack speed approaches the Rayleigh velocity (m —> l/h^), the ratio becomes zero since at that 
speed the classical ERR is infinite. It is worth noting that this trend for the ratio of the ERRs 
was also observed in the case of a mode II crack propagating with a sub-Rayleigh speed in a 
plane triangular-cell lattice (Kulakhmetova et al., 1984; Nieves et al., 2013). In fact, the plane 
triangular-cell lattice employed in these works consists of point particles (connected by massless 
elastic bonds) which interact only with forces and not internal moments, thus, corresponding 
to the case of zero micro-rotational inertia (ho = 0) in the present study. Moreover, in the 
latter studies, the global ERR for the crack propagating through the homogenized medium is 



naturally identified as the classical (far- field) ERR (Eq. (127) in our study), whereas the local 
ERR for the crack propagating through the lattice (i.e. the energy spent on fracture itself) 



corresponds to the ERR in couple-stress elasticity (Eq. ( 126 ) in our study) evaluated at the 



crack-tip (see Fig. 10). 

In the case ho > /iq, a different trend is observed for the variation of the ratio J d /J^ (Fig. 
|l2p). Specifically, as ho increases, the ratio becomes eventually a monotonically increasing 
function of the normalized speed. However, in all cases, the ratio remains below unity. More- 
over, since ho > h^, the pertinent Rayleigh speed is defined as mjj = 1/ho, which is evidently 



less than the classical Rayleigh velocity (1/ho) 
terminate before they reach the zero value. 
Finally, Figure 



Therefore, the curves plotted in Figure 



12 



13 



displays the dependence of the ratio J d / J% upon the ratio of lengths 
£/L. It is observed that the ratio J / ' J^, decreases monotonically with increasing values of£/L. 
This decrease is more significant as the propagation speed approaches the Rayleigh velocity 
rriR. This finding shows that the couple-stress theory predicts a strengthening effect since a 
reduction of the crack driving force takes place as the material microstructure becomes more 
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Figure 13: Variation of the ratio of the dynamic ERRs in couple-stress and classical elasticity with respect 
to the ratio of the characteristic lengths £/L. 

pronounced. In fact, as the propagation velocity increases, the stenghtening effect is more 
evident. Also, as l/L — > 0, the J-integral in couple-stress elasticity tends continuously to its 
counterpart in the classical theory. In the case m — > 0, the steady state solution approaches 
the static solution, confirming previous results obtained by Gourgiotis et al. (2012) for the 
stationary crack. 

7 Conclusions 

In the present work, which is closely allied in scope to Mishuris et al. (2012), we examined 
the plane-strain problem of a semi-infinite mode II crack propagating steadily with constant 
sub-Rayleigh speed in a body with microstructure governed by couple-stress elasticity. The 
case of shear loading was chosen since, in principle, couple-stress effects are predominant in 
this type of deformation (Huang et al., 1999; Gourgiotis and Georgiadis, 2007; Gourgiotis et 
al., 2012). An exact full field solution was obtained by using the Fourier transform and the 
Wiener-Hopf technique. The effect of rotational micro-inertia was also considered in our study 
since previous experience with couple-stress analyses of surface waves and anti-plane crack 
problems showed that this term is important, especially at high frequencies. It is remarked 
that this is the first study in the literature that includes micro-inertial effects in a dynamic 
plane-strain crack problem in the context of a generalized continuum theory. Our goal was to 
determine possible deviations from the predictions of classical linear elasticity when a more 
refined theory is employed to attack plane-strain crack propagation problems. By including 
pertinent characteristic material lengths, the theory of couple-stress elasticity utilized here 
accounts for effects of microstructure encompassing, thus, the analytical possibility of size 
effects, which are absent in the classical theory. 

Our results differ in several important respects from the predictions of standard linear frac- 
ture mechanics. It was shown that the microstructural parameters, such as the characteristic 
material length £ and the micro-rotational inertia h, and on the other hand the maximum 
crack speed, strongly influence the fracture process near a rapidly moving crack. Indeed, it 
was observed that depending on the values of the microstructural ratio Hq = h/l, the ac- 
tual limiting crack speed under plane strain loading conditions, can be much lower than the 
classical Rayleigh wave speed predicted by the standard elastodynamic theory. Experimental 
findings (Livne, et al., 2010) shed light to the influence of the fracture process on the limiting 
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value. Specifically, for the steady-state crack propagation problem examined here, we have 
found that both the SIF and the ERR (J-integral) depend upon the crack speed. In partic- 
ular, our results showed that as the crack speed approaches the pertinent Rayleigh velocity 
both quantities increase significantly but remain always finite. These findings are in marked 
contrast with the classical elasticity results where it is known that: (i) the respective SIF for 
the steady-state crack propagation problem does not depend upon the crack-tip speed, and 
(ii) the ERR becomes infinite as the propagating crack reaches the Rayleigh velocity. Further, 
as in the case of the stationary crack, it was observed that the ERR decreases monotonically 
with increasing values of the ratio of the characteristic material length I over the pertinent 
geometrical length L. This means that a decrease of the value of J-integral is noticed in com- 
parison with the classical theory. Therefore, the couple-stress predicts a strengthening effect 
since a reduction of the crack driving force takes place as the material microstructure becomes 
more pronounced. 

In light of the above, we conclude that couple-stress elasticity can provide a clearer picture 
of the failure process near a rapidly moving crack tip than classical elasticity. Motivated by 
experimental observations (Rosakis et al., 1999), we intend, in a future work, to extend the 
present study to the case of intersonic mode II crack propagation with the aim of gaining 
further insight into the physical mechanisms governing the fracture process of microstructured 
materials. 

Appendix A Branch cuts for the functions x{ z )i P{z) and 7(2) 

The complex function x( z ) nas f° ur branch points at ±i6i and ±^62, where (61,62) are given 



in Eq. (72). In the case ho < 1, (61,62) are always real, consequently, the branch points are 



located along the imaginary axis (Fig. Ala). On the other hand, when Iiq > 1, (61,62) are 
complex and the branch points are symmetrically located with respect to the real axis at the 
four quadrants of the complex z-plane (Fig. A. lb). Depending on the parameters /iq and 
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Figure A.l: Branch cuts for the function x(. z ) 

to, the branch cuts for the functions (3{z) and j(z) are illustrated in Figures A. 2 and A. 3, 
respectively, with e being a real number such that e —> +0. In fact, introducing e facilitates 
the introduction of the branch cuts for the function 7(2). It is noted that, in all cases, to < m,R 
with ran = min{l//io, 1/6q} in order for the crack to propagate with a sub-Rayleigh speed. 
Finally, we remark that the specific introduction of the branch cuts secures that the functions 
(x, /?,7) are always single-valued and positive along the real axis. 
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Figure A. 2: Branch cuts for the function /3(z) 
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Figure A. 3: Branch cuts for the function 7(2) 
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Appendix B Factorization of the kernel function N(z] 

The kernel N(z) is given by the expression 
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where the complex function 6 = 9{z) is defined in (85). It is noted that the kernel function 



N{z) does not have any poles or zeros in the finite complex domain in the sub-Rayleigh 
regime, only branch points (BPs) which depend upon the values of (hn,m). In particular, we 
distinguish the following three cases: 



Case I. hn < 1 and m < — : ; then the BPs are: and ±ibn. 



hn 



Case II. hn < 1 and m > 
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; then the BPs are: and ±i&i. 







Case III. hn > 1; then the BPs are: and ±ibn. 
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Note, that m < itir, always. Furthermore, for all cases (I)-(III), when Re(z) = and e < 
\z\ < I 6q I, the real and imaginary parts of N(z) take the following form 
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whereas for case II, when Re(z) = and |brj| < \z\ < \bi\, the real and imaginary parts of N(z) 
become 
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Figure B.l: Branch cuts and integration paths for the factorization of the kernel function N(z) 



Next, taking into account that N(z) — > 1 as \z\ — > oo and employing Jordan's Lemma, we 
evaluate the functions N + (z) and N~(z) defined in Eqs. (91 ) and (92), by closing the original 
integration paths C u and Cd, which extend parallel to the real axis in the complex z-plane, with 
large semi-circles at infinity on the upper and lower half-planes respectively, as it is shown in 
Fig. 4. Then, an alteration of the integration contour along with the use of Cauchy's theorem, 
allows taking as equivalent integration paths the contours C' u and C' d , around the pertinent 
branch cuts of N(z) (Fig. Bl). In particular, we obtain for the function N + (z) 



Cases I and III 



f 1 // 



In 



+ 



>C U . 



f i ( r ibo -i 

= exp< — / tan 

/—is 
tan -1 
-i&0 

-expi-i r" tan-4 ImiV(C) 
-expj J^ tan ^^ 






-ImiV(C)' 

ReiV(C) 



c- 



'lmiV(C) 
ReiV(C) 

d( ] 



d( 

c- 



z 



Case II 



N~ 



+ 



i r ibi _i 

tan 



-ibi 



—is 



+ / + 



la J- ibl 



^\ ln[N(Q} d( .} 



= eXP \-vri 
where the function tan () is defined as 



'ImJV(C)l d( 
.ReiV(C)J ( + z. 



tan — = tan - — - + < 

x \x\ 



0, for x > and y > 

7r, for x < 

27r, for x > and y < 



(B.4) 



(B.5) 



(B.6) 



31 



Also, it is noted that the contour integrals around the branch points of iV(£) in Eqs. (B.4) 
and (B.5) are all zero. Indeed, by writing £ = Cfc + r oe , where £& = ±ik are the branch points 



of N(Q, with k = (e, bo, &i), and by taking into account that (£ — £&) • In N(() — > uniformly 
as tq — > 0, it can readily shown that 
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Finally, letting e — > and combining (B.4) and (B.5), we obtain a single formula for N + (z 
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where the constant a depends on the branches of N{z) and is defined as 
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with JC( ) being the Heaviside step function. Similarly, integrating along C' u , we evaluate the 



function ./V (z) given in (92). 



Appendix C The SIF in the equilibrium case 



In the limit m — > 0, we have, according ( |73| ) and ( |B.9| ), that bo = 1 and a = 0. In addition, 
the real and imaginary parts of N(z) become, respectively 
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for < Im(z) < 1 and Re(z) = 0. In view of the above, the function N + (z) takes the following 
form 
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which after changing the variable of integration from £ to iq and by taking into account the 
properties of the inverse tangent function, yields 
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Finally, upon substituting z = il/L into the above expression, we evaluate the equilibrium SIF 



through the relation (117), which agrees exactly with the SIF given previously by Gourgiotis et 



al. (2012) for the stationary mode-II crack (see Eq. (126) in the cited work, with tq = To/L). 
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